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ABSTRACT 

While a white dwarf is, from a theoretical perspective, the most plausible primary star in Type la 
supernova (SN la), many other candidates have not been formally ruled out. Shock energy deposited in 
the envelope of any exploding primary contributes to the early SN brightness and, since this radiation 
energy is degraded by expansion after the explosion, the diflFusive luminosity depends on the initial 
primary radius. We present a new non-detection limit of the nearby SN la 201 Ife, obtained what 
appears to be just 4 hours after explosion, allowing us to directly constrain the initial primary radius, 
Rp. Coupled with the non-detection of a quiescent X-ray counterpart and the inferred synthesized 
^^Ni mass, we show that Rp < 0.02Rq (a factor of 5 smaller than previously inferred), that the average 
density of the primary must be pp > 10^ gm cm~^, and that the effective temperature must be less 
than a few xlO^ K. This rules out hydrogen burning main sequence stars and giants. Constructing 
the helium-burning main sequence and carbon-burning main sequence, we find such objects are also 
excluded. By process of elimination, we find that only degeneracy-supported compact objects — WDs 
and neutron stars — are viable as the primary star of SN2011fe. With few caveats, we also restrict the 
companion (secondary) star radius to ^ 0.1 Rq, excluding Roche-Lobe overflowing red giant and 
main-sequence companions to high significance. 

Subject headings: supernovae: general — supernovae: individual (201 Ife) — white dwarfs 



1. INTRODUCTION 

While the nature of the explosion that leads to a Type 
la supernova (SN la) — detonation (Woosley et al. 1986; 
Woosley & Weaver 1994; Livne & Arnett 1995; Fink 
et al. 2007), defiagration (Nomoto et al. 1976), or both 
(Khokhlov 1991) — and the process that leads to the ex- 
plosion trigger are not well known, it is commonly as- 
sumed that an SN la is powered by the explosion of a 
white dwarf (WD) at a pressure and temperature suf- 
ficient to ignite carbon (Nomoto 1982; Iben & Tutukov 
1984). All viable SN la models of the progenitor system 
include a companion (secondary) star that transfers mass 
(either steadily or violently) to the WD (Nomoto et al. 
1997; Podsiadlowski et al. 2008). Single-degenerate chan- 
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nels involve transfer of mass from a giant, main sequence, 
or helium star (Whelan & Iben 1973; Nomoto 1982; Mu- 
nari & Renzini 1992; van den Heuvel et al. 1992; Han & 
Podsiadlowski 2004). Double-generate channels involve 
the merger of two WDs (Webbink 1984; Iben & Tutukov 
1984). 

There is considerable circumstantial evidence that the 
primary is a C+0 WD (Hoyle & Fowler 1960), but until 
recently (Nugent et al. 2011; see also Brown et al. 2011 
for a similar analysis) there have been very few direct 
constraints. The evidence is (1) neither hydrogen nor 
helium is seen in SNe la (Leonard 2007), and few as- 
trophysical objects lack these elements, (2) the elements 
synthesized in SN la are consistent with the fusion chain 
leading from carbon going up to the iron peak, (3) de- 
generate objects can result in runaway, explosive burn- 
ing, (4) the energy gained from burning a WD matches 
that seen in a SN la, (5) simulations of the explosions of 
CO WD stars are successful at reproducing SN la spec- 
tra (e.g., Nomoto 1982). However, the discovery of a 
class of SNe la that require a WD mass above the Chan- 
drasekhar limit (Howell et al. 2006) has caused some to 
question whether a WD is involved in the explosion after 
ah (Taubenberger et al. 2011). 

A normal SN la, SN2011fe, was discovered in the 
Pinwheel Galaxy (MlOl) more than two weeks before 
it hit maximum brightness on 12 Sept 2011 UT (Nu- 
gent et al. 2011). To date, SN2011fe provides some 
of the best constraints on the progenitor system of an 
SN la: coupled with a well-measured distance modulus 
to MlOl (DM=29.05 ± 0.23 mag; Shappee & Stanek 
2011), the non-detection of a quiescent counterpart at 
optical, infrared, mid-infrared and X-ray wavebands were 
used to placed strict limits on the nature of the progen- 
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itor system. With optical imaging reaching ~100 times 
fainter than previous hmits, Li (2011), Nugent et al. 
(2011), and Horesh et al. (2011) showed that Roche- 
lobe overflow red giants as the secondary star were ex- 
cluded; He-star + WD progenitor systems were also 
largely excluded. Nugent et al. (2011) placed constraints 
on double-degenerate models based on the non-detection 
of early-time emission from shock interaction with the 
disrupted secondary WD material (Shen et al. 2011). 

Rather than focus on the progenitor system as a whole, 
in this Letter we investigate what can be gleaned about 
the primary star, the body directly responsible for pow- 
ering the SN. Nugent et al. (2011) have previously noted 
that the primary size was small (Rp < O.IRq) based on 
considerations of the early-time light curve. They con- 
cluded, based also on carbon and oxygen observed in 
the early-time spectra, that a C+0 WD was the likely 
primary. Here, coupling a new, more stringent radius 
measurement and explicitly discussing mass constraints 
(based on explosive yield) and temperature constraints 
(from quiescent non-detection) we narrow the parameter 
space for the primary even further. Whereas the Nugent 
et al. (2011) radius constraint was insufficient to rule out 
carbon-burning main sequence stars, our results appear 
to exclude such bodies. By process of elimination, we 
flnd that a WD (or a neutron star) are the only allow- 
able primary candidates. 

2. PRIMARY CONSTRAINTS 
2.1. Mass 

As a manifestly "normal" SN la (spectroscopically 
and in peak brightness), SN2011fe is expected to show 
the characteristic decline of the half-life of the radioac- 
tive process ^^Co^^^Fe, which in turn suggests at least 
0.5Mq of synthesized ^^Ni powered the early light 
curve (Hoeflich & Khokhlov 1996; Milne et al. 1999). We 
thus consider Mp^nm = O.SM© as a conservative lower 
limit to the mass of the primary; this value is below the 
lower mass limit (Mp^nm = 0.7) for sub-luminous SN la 
models (Sim et al. 20i0; Woosley & Kasen 2011). While 
a Chandrasekhar mass, Mch = IAMq is typically in- 
voked for the primary, there are no stringent upper limits 
on the primary mass. To accommodate so-called super- 
Chandrasekhar events, we thus consider a conservative 
mass range for the primary is Mp = 0.5 to 3Mq. 

2.2. Radius 

The earliest observations of SN2011fe can be used to 
constrain the radius of the primary star. The detec- 
tion at what was determined (from a remarkable several- 
day t'^ behavior of the optical light curve) to be 11 hr 
post explosion, placed a constraint on the primary to 
be Rp ^ O.IRq Nugent et al. (2011). Starting about 
7.5 hours before the PTF discovery image, we had for- 
tuitously acquired a series of images of MlOl, covering 
the position of SN2011fe for one hour; the data were 
obtained on PIRATE (Physics Innovations Robotic As- 
tronomical Telescope Explorer) on the The Open Univer- 
sity's 0.4m telescope in Mallorca (Holmes et al. 2011)^^. 

PIRATE is mainly deployed as part of the SuperWASP con- 
sortium for follow-up photometry of exoplanet transit candidates, 
and for a routine nova monitoring program of M31. 



We analyzed these images and found no significant excess 
flux at the SN location (Fig. 1). Since the filter system 
was clear we translated the non-detection to a ^-band 
magnitude equivalent under the assumption of a range 
of blackbody temperatures of the shock. For blackbody 
temperatures in the range 3000-150000 K, we find a ro- 
bust upper limit of g = 19.0 mag at 5 <j (at a mean 
time of 3.92 hr post explosion). We use this new non- 
detection to constrain the primary radius under several 
different shock models. 
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Figure 1. (top) Stacked PIRATE/Clear image obtained starting 
at ^0 + 3.38 hr (mid-point to + 3.92 hr) centered on the position 
of SN2011fe; no significant fiux at the SN position is detected, 
(bottom, left) Zoomed images near the SN position, showing, from 
left to right, a fake star with g = 18.5, 19.0 and 19.5 mag. The 
star was created using a postage stamp image of an isolated bright 
star on the stacked image with known SDSS ^f-band magnitude 
and a color consistent with Teff = 4000 K. A shock with g = 19.5 
mag would have been marginally detected. The compact source 
labelled "A" is SDSS J140306.16+541706.5, with g = 18.20 and 
g — r = 1.16. (bottom, right) Same as at left, but using a known 
SDSS star with Tgff = 7000 K. A shock with g = 19.0 mag would 
have been marginally detected. 



Initially, the explosion shock wave deposits radiation 
energy throughout the stellar envelope, which subse- 
quently diffuses out and contributes to the supernova 
brightness. Because the radiation energy is lost to adia- 
batic expansion after the explosion, the luminosity from 
shock heating depends on the initial stellar radius. A 
number of analytical models for the early luminosity have 
been constructed under the assumption of spherical sym- 
metry (e.g.. Chevalier 1992; Chevalier & Fransson 2008; 
Piro et al. 2010; Kasen 2010; Nakar & Sari 2010; Rabi- 
nak & Waxman 2011; Rabinak et al. 2011). The studies 
differ in their assumptions of e.g., the initial ejecta den- 
sity and pressure profiles, the nature of the opacity, and 
in the treatment of radiative diffusion. Nevertheless, the 
predicted early light curves tend to be quite similar, and 
have been found to agree with numerical calculations to 
within a factor of ~ 2 (Rabinak et al. 2011; Kasen 2010). 
For example, Rabinak et al. (2011) find analytic expres- 



A Compact Degenerate Primary-Star Progenitor of SN2011fe 



3 



sions for the luminosity and effective temperature: 



L{t) = 1.2 X 10' 
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p 77^0.85 
^10^51 



t~^-^^ ergs s"^ 



Teff(t) =4105 



p 

^0.2 
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f0.022 
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where Rio is the progenitor radius Rp/ 10^^ cm, £^51 is 
the explosion energy ergs, Mc is the mass in units 

of Mch, and ^^0.2 is the opacity Hi/0.2 cm^g~^. These ex- 
pressions assume a constant opacity, which is appropriate 
when electron scattering dominates in fully ionized C/0 
eject a. The form factor /p depends on the density pro- 
file of the primary star, and has estimated values in the 
range 0.031 and 0.13 (Calzavara & Matzner 2004; Rabi- 
nak et al. 2011). 

We constructed theoretical light curves of the early op- 
tical luminosity assuming the spectrum was given by a 
blackbody with Teff, and the flux at the effective wave- 
length of the filter (given the temperature) was consis- 
tent with the non-detection. Figure 2 shows the results 
for a selection of different analytical models, assuming 
E^i/Mc = 1 (constant explosive yield per unit mass), 
fp = 0.05, /^o.2 = 1 and a variety of values of Rp. The 
PIRATE observation at 4 hours is the most constraining 
data point, which limits Rp < 0.02Rq . Table 1 sum- 
marizes the detailed radius constraints under different 
assumptions of progenitor mass and under the different 
models. 

The expressions we have used for the early luminosity 
hold only under the assumption that radiation energy 
dominates in the post-shock ejecta. In fact, the diffusion 
wave will eventually recede into higher density regions of 
ejecta where gas pressure dominates. The luminosity is 
then expected to drop suddenly; Rabinak et al. (2011) 
show that, for constant opacity, the time of this drop is 
proportional to Rp, which effectively limits the minimal 
progenitor radius that we are capable of probing. From 
their expression for tdrop we find this minimal radius to 
be 



Rr. 



0.013 t4hi?5"i° ''^c 'Vo.ol R' 



(2) 



where = ^4 hours. The value of i^min is just smaller 
than our limits on Rp determined in Table 1, suggesting 
that the breakdown of radiation energy domination is not 
likely to undermine our results. 

The early photometry of SN2011fe also tightly con- 
strain the nature of a possible companion star. The in- 
teraction of the supernova ejecta with a companion star 
produces emission which depends linearly on the sep- 
aration distance (Kasen 2010). This emission will be 
anisotropic and vary by a factor of ~ 10 depending on 
the orientation. Assuming the companion star in Roche 
Lobe overflow, such that its radius is < 1/2 of the sep- 
aration distance, and that the observer's viewing angle 
is unfavorable (such that the light curve is fainter by a 
factor of 10 from its maximum) our data restricts the 
companion star radius to i?c ^ 0-1 ^o- Unless the time 
since explosion for the PIRATE data is vastly underesti- 
mated (by ^ day), this apparently excludes Roche-Lobe 
overflowing red giant and main sequence companions to 
high significance. 

Temperature- Radius: Non-detections of a quiescent 
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Figure 2. Absolute g-band magnitude versus time since explosion 
in three theoretical models for the early time evolution of Type la 
SNe. Shown is 4 hr, 5 cr non-detection discussed in §2.2 and the 
first two detections from PTF (Nugent et al. 2011). The black line 
shows the L oc behavior of the early-time light curve detected 
in PTF, consistent with the non-detection. For the Kasen (2010) 
companion interaction model, R denotes the separation distance 
between the two stars, and the light curve is shown for an ob- 
server aligned with the collision axis, which produces the brightest 
observed luminosity. 



counterpart in Hubble Space Telescope (HST) imaging 
yield a specific luminosity (L^) constraint at certain op- 
tical frequencies. With the assumption of a spectrum 
of the primary, these limits can be turned into a limit 
on the bolometric luminosity (L). Li (2011) considered 
mostly spectra of an unseen secondary, using model input 
spectra of red giants to derive L constraints. For a high 
effective temperature primary, here we consider a black- 
body as the input spectrum and solve for the bolometric 
luminosity and effective radius using Stephan-Boltzmann 
law (see also Liu et al. 2011 for a similar analysis). We 
perform a similar analysis with the Chandra X-ray non- 
detection, convolving different input blackbody spectra 
to find a radius limit. At 10^ K, for example, the limits 
(1,2, and 3 a) are 1.2 x 10"^ Rq, 1.5 x 10"^ Rq, and 
1.8 X 10-^ Rq. 

In Figure 3, we show these primary- star constraints as 
a function of effective temperature and radius. Primary 
stars with average density less than pp = 10^ gm cm~^ 
and effective temperatures larger than 10^ K (at pp = 
10^^ gm cm~^ ) are excluded. 

3. COMPARISONS TO PRIMARY CANDIDATES 

Accepting O.SM© as a conservative lower limit for 
the primary mass, low-mass main-sequence stars, brown 
dwarfs, and planets are not viable. In Figure 3, we 
show the main sequence of stably H-burning stars with 
mass 0.5, 1, 1.4, and 3 Mq. The hydrogen main se- 
quence, shown using solar- metallicity isochrones from 
Marigo et al. (2008), is excluded as the SN2011fe pri- 
mary. Accepting the radius constraints, giants (not plot- 
ted) are also excluded for the primary of SN2011fe. 

We also constructed idealized He- and C-burning main 
sequence stars with the stellar evolution code MESA 
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Figure 3. Constraints on mass, effective temperature, radius and 
average density of the primary star of SN2011fe. The shaded red 
region is excluded from non-detection of an optical quiescent coun- 
terpart in Hubble Space Telescope (HST) imaging. The shaded 
green region is excluded from considerations of the non-detection 
of a shock breakout at early times, taking the least constraining Rp 
of the three model in Table 1. The blue region is excludes by the 
non detection of a quiescent counterpart in Chandra X-ray imag- 
ing. The location of the H, He, and C main sequence is shown, 
with the symbol size scaled for different primary masses. Several 
observed WD and NSs are shown. The primary radius in units of 
Rq is shown for Mp = IAMq. 

(Paxton et al. 2011).^^ Uniform composition stars 
(-^4He = 0.98 and solar metallicity for the He stars, and 
Xi2c = Xi60 = 0-5 for the C stars) were relaxed un- 
til they reached stable non-degenerate equilibrium con- 
figurations. Up-to-date neutrino loss rates (Itoh et al. 
1996), opacities (Iglesias & Rogers 1996), equation of 
state (Timmes & Swesty 2000; Rogers & Nayfonov 2002), 
and nuclear reactions (Caughlan & Fowler 1988; Angulo 
et al. 1999) were employed, including recently updated 
triple-a (Fynbo et al. 2005), a+^^C (Kunz et al. 2002), 
and ^^C+^^C (Casques et al. 2005) rates. These results 
are plotted in Figure 3 for He star masses of 0.5, 1.0, 1.4, 
and 3.0 Mq, and C star masses of 1.0, 1.4, and 3.0 M©; 
stable configurations of C stars supported by nuclear 
burning do not exist below 1.0 Mq (Boozer et al. 1973). 
Comparisons to previous He star calculations (e.g.. Di- 
vine 1965; L'Ecuyer 1966) and C star calculations (e.g., 
Sugimoto et al. 1968) match to ^ 10% in the radius and 
effective temperature. Evidently, He-main sequence stars 
and C-burning main sequence stars are also excluded as 
the SN2011fe primary. 

Several WDs in eclipsing binary systems, with mea- 
sured temperature, radii, and mass, are shown in Fig- 
ure 3 (SDSS1210: Pyrzas et al. 2011; V471 Tau: Pyrzas 
et al. 2009; Sirius B: Barstow et al. 2005). We also de- 
pict isolated neutron star RX J185635-3754 (M ^ Mq), 
which has strong observed constraints on temperature 
and density (Pons et al. 2002). These systems are all 
allowed by our constraints. Note that SDSS 1210 (M = 
O.415±O.1M0; i? = O.O159±O.OO2i?0; 1.3-1.6x10^ 

http://mesa.sourceforge.net/ (version 3611) 



Table 1 

Primary Radius and Density Constraints 



Mp 


D a 
-fi'P,max 


avg. pp 


-^shock,max 


^shock,max 


Mq 


Rq 


gm cm~^ 


erg s"-*- 


K 




Shock breakout — 


Rabinak, Livne, &; Waxman (2011)^ 


0.5... 


0.022 


6.29 X 10^ 


4.50 X 10^9 


6024 


1.0... 


0.020 


1.77 X 10^ 


4.48 X 10^9 


6043 


1.4... 


0.019 


2.93 X 10^ 


4.47 X 10^9 


6053 


3.0... 


0.017 


9.16 X 10^ 


4.46 X 10^9 


6075 




Ejecta Heating Secondary 


— Kasen (2010)^ 


0.5... 


0.038 


1.26 X 10^ 


3.94 X 10^9 


8048 


1.0... 


0.027 


7.05 X 10^ 


3.96 X 10^9 


7388 


1.4... 


0.023 


1.58 X 10^ 


4.00 X 10^9 


7103 


3.0... 


0.017 


9.29 X 10^ 


4.18 X 10^9 


6531 


Shock breakout — Piro et al. (2010)^ 


0.5... 


0.016 


1.73 X 10^ 


5.27 X 10^9 


12110 


1.0... 


0.019 


2.07 X 10^ 


5.26 X 10^9 


12091 


1.4... 


0.021 


2.26 X 10^ 


5.26 X 10^9 


12082 


3.0... 


0.025 


2.75 X 10^ 


5.25 X 10^9 


12062 



^ 5 cr limit assuming the 4 hr non-detection (see text) and shock 

opacity k = 0.2 cm^ g~^- 

^ Assumes fp = 0.05, E^i/M^ = 1 

^ The radius derived is the separation distance and the limit de- 
rived is assuming the brightest possible viewing angle. The ra- 
dius limit comes from the requirement that primary size must be 
smaller than the semi-major axis of the binary. 
^ Using their eqns. 35 & 36 but corrected by a factor of 7~'^^^ 
(L) and 7~^/^ (Tgff) to fix the improper scalings. 

gm cm~^) is a He WD and has a mass less than our 
suggested Mp > 0.5 constraint. It is also marginally ex- 
cluded by the Piro et al. (2010) model. 

4. SUMMARY AND DISCUSSION 

We have placed limits on the average density, effec- 
tive temperature, and radius of the primary (explod- 
ing) star of the Type la SN SN2011fe. We consider the 
g = 19.0 mag non-detection as conservative and, assum- 
ing the t'^ behavior accounted for some of the flux at 4 
hours, the flux from the shock inferred from this non- 
detection would necessarily have been even smaller than 
that derived from g = 19.0 mag. In this respect, we 
take the radius constraint of the primary to be very con- 
servative with the important proviso: if the explosion 
time was significantly earlier than the time inferred from 
the fit (Nugent et al. 2011), the radius constraints 
are less stringent. In particular, if the PIRATE obser- 
vations occurred at to + 28 hr instead of to + 4 hr, then 
Rp ^ 0.2Rq (still sufficient to rule out H and He main- 
sequences but not the C MS). Another important caveat 
is that the radius constraint requires shock heating, natu- 
rally expected with a deflagration-detonation transition 
(Khokhlov 1991)^^. A pure deflagration of a WD that 
does not produce a strong shock would not exhibit the 
early-time behavior of the models presented in Table 1; 
however, pure deflagration is disfavored on nucleosyn- 
thetic grounds (Nomoto et al. 1984; Thielemann et al. 
1986). With these caveats aside, this is the most strin- 

"Double-detonation" scenarios (Woosley Sz Weaver 1994) 
could also lead to heating, but with (slightly) different heating 
than considered by the models in §3. 
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gent limits on the primary radius and temperature of an 
SN la reported to date. 

Clearly the density and temperature in the core of a 
primary are higher than the reported constrained quan- 
tities; and since at high density and high central temper- 
ature a star may be supported by pressure other than 
that associated with fermionic degeneracy pressure, we 
cannot formally exclude all non-degenerate stars. How- 
ever, by process of elimination we find that only compact 
degenerate objects (WD, NS) are allowed as the explo- 
sive primary. This statement comes from considerations 
almost orthogonal to the traditional spectral modeling in 
Type la SNe that are invoked to claim WDs as the ex- 
ploding primary Since the explosive nucleosynthetic 
yield from phase transition of an NS to a quark-star is ex- 
pected to be very small (< O.IMq in r-process elements) 
and is unlikely to produce light elements (Jaikumar et al. 
2007), a NS primary is disfavored. 

We thank Eliot Quataert, Alex Filippenko, Lars Bil- 
sten, Weidong Li, William Lee, and Philipp Podsiad- 
lowski for helpful discussions. 
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■•^^ As noted in Nugent et al. (2011), the early observations of C 
and O in the spectrum are highly suggestive that the primary was 



a C+O WD. 



